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Abstract 

We analyzed the experimental data of the pre-scission neutron multiplicity in con- 
nection with fission fragments in the reaction 58 Ni+ 208 Pb at the incident energy 
corresponding to the excitation energy of compound nucleus £'*=185.9 MeV, which 
was performed by DeMoN group. The relation between the pre-scission neutron 
multiplicity and each reaction process having different reaction time is investigated. 
In order to study the fusion-fission process accompanied by neutron emission, the 
fluctuation-dissipation model combined with a statistical model is employed. It is 
found that the fusion-fission process and the quasi-fission process are clearly distin- 
guished in correlation with the pre-scission neutron multiplicity. 

Key words: superheavy elements, fluctuation-dissipation dynamics, fusion-fission 
process, quasi-fission process, pre-scission neutron multiplicity 



1 Introduction 

Many experiments on the induced fission at near and below Coulomb barrier 
in superheavy-mass region have been done and the mass and kinetic energy 
distributions of fission fragments were measured (Q). In these phenomena, 
the reaction process is classified into several characteristic processes, they are 
the fusion-fission process (FF), the quasi-fission process (QF) and the deep 
inelastic collision process (DIC). Depending on the shell effect of the compound 
nucleus or the composite nucleus, the deep quasi-fission process (DQF) can 
be seen in a certain reaction system (0). In the DQF process, mass symmetric 
fission fragments are observed but no compound nucleus is formed. 
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Each process has its own characteristic reaction time from the contact of 
the colliding partner to the scission point. This was shown by the dynamical 
calculation for the time development of the nuclear shape in terms of the 
Langevin equation (Q). The different reaction time means that each reaction 
process associates the different pre-scission neutron multiplicity. 

We undertake to extend our model discussed in reference (|2j), by taking into 
account the effect of neutron emission. We combine the Langevin calculation 
with the statistical model which calculates the neutron emission. We apply our 
model to the recent experiment, in which the pre-scission neutron multiplicity 
correlated with the mass distribution of fission fragments has been measured 
in the reaction 58 Ni+ 208 Pb at the incident energy corresponding to the exci- 
tation energy of compound nucleus E* = 185.9 MeV. This experiment was 
done by DeMoN group (J3j). The feature of distribution of pre-scission neutron 
multiplicity v n has two distinct components, one located around v n = 4 and 
the other around v n = 8. This structure may be sign of two mechanisms hav- 
ing two different life-time of the composite system, they are, the QF process 
and the FF process. By our model calculation, it can be clearly shown the two 
components of the neutron multiplicity are corresponding to the above two 
processes. 

In this paper, we show theoretically the usefulness of the methods of classifying 
the dynamical process on the basis of the pre-scission neutron multiplicity 
correlated with the mass distribution of fission fragments. This approach is 
a powerful tool for investigating fusion-fission mechanism and for estimating 
precise cross section relevant to the problem on superheavy elements synthesis. 

In section 2, we briefly explain our framework for the study and the model. 
We discuss the pre-scission neutron multiplicity in correlation with fission 
fragments in section 3. We also present how the pre-scission neutron multi- 
plicity depends on the strength of the friction, the level density parameter 
and the neutron binding energy used in the statistical model. In section 4, we 
present a summary and further discussion to clarify the reaction mechanism 
in superheavy mass region. 



2 Model 



Using the same procedure as described in reference (|2j) to investigate the 
fusion-fission process dynamically, we use the fluctuation-dissipation model 
and employ the Langevin equation. We adopt the three-dimensional nuclear 
deformation space given by two-center parameterization (0; 0) . The neck pa- 
rameter e is fixed to be 1.0 in the present calculation, so as to retain the 
contact-like configuration more realistically for two-nucleus collision. The three 
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collective parameters involved in the Langevin equation are as follows: zq (dis- 
tance between two potential centers), 5 (deformation of fragments) and a 
(mass asymmetry of the colliding nuclei); a = (Ax — A 2 )/ (A\ + A 2 ), where Ai 
and A 2 denote the mass numbers of the target and the projectile, respectively. 

The multidimensional Langevin equation is given as 



dq 



,lt ( m ~% P » 



(m x ) pjp k - 7^ (m 1 ), k Pk + 9ijRj(t), 



dpi dV 1 d 

~dt = ~dq~~2dq~ i y'" h^' J " K rij y hk 

where a summation over repeated indices is assumed. denotes the deforma- 
tion coordinate specified by z , 5 and a. Pi is the conjugate momentum of 
V is the potential energy, and and 7^ are the shape- dependent collective 
inertia parameter and dissipation tensor, respectively. A hydrodynamical in- 
ertia tensor is adopted in the Werner- Wheeler approximation for the velocity 
field, and the wall-and-window one-body dissipation is adopted for the dissi- 
pation tensor 0; 0). The normalized random force Ri(t) is assumed to be 
white noise, i.e., (Ri(t))=0 and (Ri(tx)Rj(t2)) = 2SijS(ti — t 2 ). The strength 
of random force is given by y^T = J2k9ij9jk, where T is the temperature 
of the compound nucleus calculated from the intrinsic energy of the compos- 
ite system as E int = aT 2 , with a denoting the level density parameter. The 
temperature-dependent potential energy is defined as 

V(q,l,T) = V DM (q) + + V shell (q)$(T), (2) 

V DM (q) = E s (q) + E c (q), (3) 



where I(q) is the moment of inertia of a rigid body at deformation q, V s h e ii is 
the shell correction energy at T = 0, and Vdm is the potential energy of the 
finite-range liquid drop model. Eg and Eq denote a generalized surface energy 
(0) and Coulomb energy, respectively. The centrifugal energy arising from the 
angular momentum I of the rigid body is also considered. The temperature- 
dependent factor $ is parameterized as <3> =exp{— aT 2 / E^} following the work 
of Ignatyuk et al. (ficfr). The shell dumping energy E d is chosen to be 20 MeV. 
The intrinsic energy of the composite system E int is calculated for each tra- 
jectory as 

E int = E *~\ - V(q,l,T), (4) 



where E* denotes the excitation energy of the compound nucleus, and is given 
by E* = E cm — Q with Q and E cm denoting the Q— value of the reaction and 
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the incident energy in the center-of-mass frame, respectively. 



We take into account neutron emission in the Langevin calculation during the 
fusion-fission process. The present work is the first theoretical attempt to in- 
vestigate neutron evaporation in the dynamical calculation of the fusion-fission 
process in the superheavy-mass region. Here, we assume that no neutron is 
emitted in the approaching process. Using the procedure outlined by Frobrich 



et al. in which the Langevin calculation and the statistical model were 
combined, the emission of neutrons has been coupled to the three-dimensional 
Langevin equation for the fusion-fission process. Due to the strong friction in 
this model, we assume that the kinetic energy of relative motion quickly dis- 
sipates into the intrinsic energy after contact. The thermal equilibrium in the 
composite system in this calculation is realized immediately within the time 
order of 10~ 22 s (Jlih . We also assume that the particle emissions in the com- 
posite system are limited to neutron evaporation in the neutron-rich heavy 
nuclei. 

The width of neutron emission is expressed as 



2tt Pa (E, 



int ) 



f J2(2l + l)Ti{e) PA ^{E int -B n - e)de, (5) 
o « 



where g n is the spin multiplicity of the neutron. B n and e denote the neutron 
separation energy and the kinetic energy of the emitted neutron, respectively. 
The transmission coefficient is presented by TJ. The nuclear level density is 
denoted by p(E) with suffixes A and A — 1 for the compound nucleus and 
the residual one, respectively. R is the radius of the compound nucleus. Using 
the sharp cut off treatment for the angular momentum distribution of 7], we 
express Eq. (5) as follow: 

1 _'///,'- i; ""'~ 7J " 



r « = 9n- 77; — r— T2— / e PA-i(E int - B n - e)de. (6) 



where p is the neutron mass (0 EH). In the case of the excitation energy 
corresponding to a high incident energy compared with the Bass barrier height 
(jlih . as in the present case, the integration can be performed approximately, 
if we assume a simple form of p{E) oc e 2y/anE for the level density flih. as 
follows: 



.— ^{E int - B n )exp{-J-?-B n }, (7) 
Tin a n U rj int 



where a n is the level density parameter. 
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The Langevin calculation is solved by using constant time steps At, and during 
each time step, the neutron emission is calculated by the following Monte Carlo 
procedure. We obtain the lifetime of neutron emission r n = h/T n from Eq.(7). 
The probability of neutron emission during the time interval At is described by 
— . If the ratio — (in the actual calculation, At < r n ) is larger than a random 
number £ (0 < £ < 1), the emission of a neutron is registered. The kinetic 
energy of the emitted neutron eW is estimated stochastically by assuming the 
neutron spectrum from the composite system with the temperature T to be 
e 1//2 exp (— e/T). The upper index (i) indicates the i-th neutron emission stage. 
Therefore, the energy loss in the neutron evaporation process is estimated as 
_ e (») —B^}. Then we set a new intrinsic energy E^ t = Ef nt ^ —e®—B$, where 
for i — 1, E^t means the initial value of Ei nt . The potential energy surface also 
changes as described by Eq.(2) as the temperature decreases. After the initial 
energy is reset, we proceed to the next time step of the Langevin calculation. 



3 Results and discussion 

Using our model, we analyze the experimental data in the reaction 58 Ni+ 208 Pb 
at i?*=185.9 MeV (0). Figure 1 shows the experimental result, which is the dis- 
tribution of the pre-scission neutron multiplicity v n in correlation with fission 
fragments whose mass number is greater than ^ — 30 and less than ^ + 30. 

An interesting feature of the distribution of pre-scission neutron multiplicity 
is its shape having two components, which are located around v n = 4 and 
v n = 8. This structure may be the sign of a simultaneous coexistence of two 
mechanisms corresponding to different life time of the composite system but 
to the phenomena giving nearly the same mass fragment in the fission process. 
The first one, defined as the QF process, would lead to the emission of about 

4 neutrons only and for the second one, associated here with the FF process 
via a compound nucleus, would be found around 8 neutron emission. In this 
case, the fission process would tail long enough time to allow the emission 
of nearly 8 neutrons. By our model calculation, we try to confirm that the 
different dynamical processes, i.e., the FF process and the QF process, are 
giving different pre-scission neutron multiplicities in spite of the phenomenon 
associated with the similar mass fragment. 

As mentioned in the previous section, in the calculation of neutron emission, 
there exist two important parameters, the level density parameter a n and the 
neutron binding energy B n . Both of these parameters depend on the defor- 
mation of the nucleus. In the dynamical calculation itself, another unknown 
factor is the dissipation energy Edi SS , which is converted from the relative 
kinetic energy of the colliding system to the intrinsic energy during the ap- 
proaching process. In this section, first we present the structure appeared in 
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Fig. 1. Experimental pre-scission neutron multiplicity associated with fission frag- 
ment measurements in the reaction 58 Ni+ 208 Pb at E*=185.9 MeV Q). 

our model calculation of the pre-scission neutron multiplicity. Secondly, we 
discuss how the results depend on the parameters mentioned above. 

3. 1 Neutron emission during fusion-fission process 

As discussed in reference (0), in order to simplify the calculation, we assume 
again that the kinetic energy does not dissipate during the approaching pro- 
cess, that is, E diss = MeV. In subsection 3.2.2 below, the sensitivity of E diss 
is discussed. In the three-dimensional Langevin calculation, we prepare 1,000 
trajectories, and at t = 0, each trajectory starts from the point of contact, 
which is defined as Rtouch = R P + Rt, where R p and R t are the radii of the 
projectile and the target, respectively. 

The level density parameter a n for a certain nuclear shape has been estimated 
by classical methods ffl6h . Fundamentally , we need a n for various deformations 
of the nucleus. According to reference (|16l). the value of a n increases with 
increasing deformation of the nucleus. In order to simplify the calculation, we 
use two different values of a n , that is, the value inside the fusion box defined 
in the deformation space as {z < 0.5,5 < 0.3, \a\ < 0.3}, and that outside 
the region, which are denoted by a^ h and a„ * , respectively. The excitation 
energy dependence of a n is not considered. The neutron binding energy B n 
also depends on the deformation, but here we use the constant value given 
for the ground state, because the fluctuation of B n in the entire deformation 
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space relevant to the present calculation is ±0.3MeV. The influence of the 
fluctuation of B n is discussed further in subsection 3.2.4. We chose the values 
of a sph and so as to reproduce the features of the experimental data, which 
are a s n ph = A/12.0 and < e/ = A/10.0. 
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Fig. 2. Pre-scission neutron multiplicity given by our theoretical calculation in the 
reaction 58 Ni+ 208 Pb at £*=185.9 MeV. The neutron multiplicity from the QF pro- 
cess and the FF process are denoted by the blue and red line, respectively. The black 
line denotes the total processes, (a) in the case of a sph = A/12.Q and afff = A/10.0 
(b) in the case of a s n ph = A/14.0 and a^ f = A/10.0. 

For the QF process, we plot the pre-scission neutron multiplicity in correlation 
with fission fragments with mass numbers greater than 4 — 30 and less than 
~ + 30. To classify the trajectory of the FF path, the definition of the fusion 
area in the deformation space is very important. Here, we define the fusion area 
(fusion box) as the inside of the fission saddle point in the system described 
above. The idea behind the definition of the fusion box is the same as that in 
reference (0). The FF trajectory is identified as that which enters the fusion 
box. 

The distribution of the pre-scission neutron multiplicity v n calculated by our 
model is shown in Fig. 2(a) with the level density parameters a sph = v4/12.0 
and = A/10.0. The pre-scission neutron multiplicity for the QF and the 
FF processes are denoted by the blue and the red lines, respectively. The black 
line shows the total multiplicity of each process. As another example, we show 
in Fig. 2(b) the result using the different level density parameter set, which 
are a s n ph = A/14.0 and = A/10.0. 
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Fig. 3. Sample trajectories projected onto z — a (5 = 0) plane at E* = 185.9 MeV 
in the reaction 58 Ni+ 208 Pb. The trajectories of the QF and the FF processes are 
denoted by blue and red lines, respectively. The potential energy surface is presented 
by the liquid drop model in nuclear deformation space for 266 Ds. The arrow denotes 
the injection point of the reaction. 

We can see clearly the two components which come from the QF process and 
the FF process. It is shown that for the large neutron multiplicity it originates 
from the FF path, and on the other hand for small neutron multiplicity it 
comes from the QF path. This means that the pre-scission neutron multiplicity 
has a strong correlation with dynamical paths. 

On the pre-scission neutron multiplicity, the odd-even oscillations appear 
clearly, due to the odd-even effect on neutron binding energies. In the system, 
the initial number of neutrons is even, so that the probability of emitting an 
even number of neutrons is larger. The calculations show the similar structure 
observed in the experimental measurements in Fig. 1. 

Next, we discuss the details on this calculation. Figures 3 and 4 show the 
potential energy surface of the liquid drop model for 266 Ds on the z — a (5 = 
0) plane and z — 6 (a = 0) plane, respectively, in the case of I = 0. This 
potential energy surface is calculated using the two-center shell model code 
ffl7l : llif). The contour lines of the potential energy surface are drawn at steps 
of 2 MeV. Similar to the deformation parameters of nuclear shape described 
in reference (0), z is defined as z = z /(R C nB), where Rcn denotes the 
radius of the spherical compound nucleus. The parameter B is defined as 
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Fig. 4. The same samples of the trajectory as in Fig. 3, but which are projected 
onto z — 5 (a = 0) plane. 

B = (3 + 5)/(3 -25). 

In Figs. 3 and 4, the position &tz = a = 5 = corresponds to a spherical 
compound nucleus. The injection point of this system is indicated by the arrow. 
The top of the arrow corresponds to the point of contact in the system. We 
start the calculation of the three-dimensional Langevin equation at the point 
of contact, which is located at z — 1.575,5 = 0.0, a = 0.564. All trajectories 
start at this point with momentum in the initial channel. That is to say, the 
initial velocity is directed in only the z direction, and the components of the 
initial velocities along the S and a directions are both assumed to be zero. 
The sample trajectories of the QF process and the FF process are shown in 
Figs. 3 and 4. The trajectories are projected onto the z — a plane (5 = 0) in 
Fig. 3 and z — 5 plane (a = 0) in Fig. 4. The trajectory of the QF process and 
the FF process are denoted by blue line and red line, respectively. 

We define the travelling time t tra v as a time duration during which the trajec- 
tory moves from the point of contact to the scission point. Figure 5 shows the 
distribution of travelling time t trav . The t trav from the QF and FF processes 
are denoted by the solid and dashed lines, respectively. On the FF process, as 
we can see in Figs. 3 and 4, the trajectory is trapped in the pocket around the 
spherical region. The trajectory spends a relatively long time in the pocket and 
it has a large chance to emit neutrons. In average, the time duration spending 
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Fig. 5. The distribution of travelling time tf rav in the reaction 58 Ni+ 208 Pb at 
E* = 185.9 MeV. The t trav from the QF and FF processes are denoted by the 
solid and dashed lines, respectively. 

in the pocket is 7 x lCT 20 sec. On the other hand, the trajectory of the QF 
process reaches to the scission point quickly. So, it has not so much chance to 
emit neutrons. The mean travelling time for the QF process is approximately 
2 x 10- 20 sec. 

Actually the distribution of t trav for each trajectory is clearly distinguished 
by the two processes as shown in Fig. 5. The time scale of the FF process is 
about 3 or 4 times longer than that of the QF process. We can see in Figs. 3 
and 4 that the length of QF path is shorter than that of the FF path. On the 
other hand, the FF trajectory is trapped at the potential pocket near spherical 
region, so that the traveling time becomes longer. 

It is interesting to investigate the correlation between the travelling time and 
a number of events for neutron multiplicity. Figure 6 shows the relation of the 
travelling time t trav and pre-scission neutron multiplicity v n . The blue dots and 
red dots are denoted the events for the QF and the FF processes, respectively. 
As shown in Fig. 6, the longer the travelling time is, the larger the pre-scission 
neutron multiplicity becomes. Also in the figure, the correlation between the 
mass distribution of fission fragments and the pre-scission neutron multiplicity 
is shown. The events from the QF process are distributing around u n = 2n, 3n 
and An, and in the time axis they are distributing around t tra v ~ 1-5 x 10 -20 
sec. In the case of the FF process, due to longer ttrav, the events are scattered 
around v n = 6n, 7n and 8n. Since the mass distribution of fission fragments 
shows symmetric in this reaction system due to high incident energy, we can 
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Fig. 6. The relation of the travelling time tt ra vi the mass distribution of fission frag- 
ments and pre-scission neutron multiplicity v n . The blue and red dots are denoted 
the events from the QF and FF processes, respectively. 

not clearly see the correlation between the mass fragment and the neutron 
multiplicity. 

3.2 Parameter dependence of the model 

Our model has several parameters which are not confirmed theoretically at 
present. We discuss the sensitivity of the parameters of our model and show 
clearly the role of each parameter. As the essential factors that influence the 
neutron multiplicity, two categories are considered. They are concerned with 
the travelling time of the trajectory and with the rate of neutron emission. 
The former comprises the strength of the friction force and the dissipation 
energy Edi SS at the point of contact. The latter comprises the level density 
parameter a n and the neutron binding energy B n . 

3. 2. 1 Friction force 7 

In the dynamical calculation, the friction force plays a very important role. The 
friction force influences mainly the travelling time of the trajectory. Moreover, 
the rate of dissipation at which the relative kinetic energy is converted into the 
intrinsic energy and the strength of fluctuation in the trajectory are controlled 
by the magnitude of the friction force. In order to elucidate the sensitivity 
of the one-body friction tensor to the trajectory calculation, the pre-scission 
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Fig. 7. The neutron multiplicity calculated using one-body friction tensor multiplied 
by factors 0.75, 1.0 and 1.5, as denoted by the dashed, solid and dashed-dotted lines, 
respectively. The neutron multiplicity of the QF and FF processes is dented by the 
black and gray lines, respectively. 

neutron multiplicity is shown in Fig. 7 for the modified strength of the one- 
body friction tensor. The multiplied factors are 0.75, 1.0 and 1.5, and the 
corresponding results are shown by the dashed, solid and dashed-dotted lines, 
respectively. The neutron multiplicities from the QF and FF processes are 
dented by the black and gray lines, respectively. 

For large friction with the factor of 1.5, the trajectory speed is slower than that 
in the case of friction with the factor of 1.0. Due to longer t tra v, the neutron 
multiplicity v n becomes larger. As a result, the peaks of neutron multiplicity 
from the QF process and the FF process shift to larger multiplicity (to the right 
in Fig. 7). On the other hand, in the case of the factor of 0.75 (smaller friction), 
the trajectory speed is faster and the neutron multiplicity distribution shifts 
to smaller v„. 



3.2.2 Dissipation energy E diss 

As the initial condition of the Langevin calculation, the intrinsic energy and 
relative kinetic energy of nuclei at the point of contact are very important. 
However, in the superheavy-mass region, it is an open question as to how 
much of the relative kinetic energy of the colliding system dissipates into in- 
trinsic energy during the approaching process. Although we assume no energy 
dissipation E diss at the point of contact in our model calculation, in order to 
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Fig. 8. The neutron multiplicity in the cases of Ediss = and 32.4 MeV, as denoted 
by solid and dashed lines, respectively. The black and gray lines denote the neutron 
multiplicity from the QF and FF processes, respectively. 



elucidate the sensitivity of Edi SS , we deal with it at that point as a parameter. 
We have discussed this subject on the basis of an analysis of the mass distri- 
bution of fission fragments in reference (jioh . Here, we discuss how the energy 
loss E diss affects the neutron multiplicity. 

Since the initial momentum in the — z direction strongly affects the trajec- 
tory, E diss is expected to play a substantial role. Figure 8 shows the neutron 
multiplicity in the cases of Edi SS = and 32.4 MeV (20% of kinetic energy dis- 
sipates), which are denoted by solid and dashed lines, respectively. The black 
and gray lines denote the neutron multiplicity from the QF and FF processes, 
respectively. 

In the case of Edi SS =32A MeV, due to the smaller initial momentum in the — z 
direction, the number of trajectories that reach the spherical area decreases. 
However, due to the odd-even effect on the neutron binding energies, such a 
variance of initial momentum is too small to significantly change the neutron 
multiplicity distribution. 

Quantitatively, the mean value of v n for the QF process increases from 3.53 
to 3.60 when 20% of the kinetic energy is converted to the internal one during 
the approaching process. On the contrary, in the same situation, the mean 
value of v n for the FF process decreases from 7.94 to 7.87. 



13 




Fig. 9. The neutron multiplicity calculated using a n = A/12.0 in the entire defor- 
mation space, denoted by the dashed lines. The solid lines denote the calculation 
using a^ h = A/12.0 and a*^ = A/10.0. The black and gray lines denote the QF 
and FF processes, respectively. 



3.2.3 Level density parameter a n 

In the statistical model calculation, the level density parameter a n is very 
important. The lifetime of neutron emission is strongly affected by a n , because 
a n appears in the exponent in the formula of the decay width, as shown by 
Eq. (7). 

First, we use the constant value a n in the entire deformation space. The value 
of a n is chosen so as to reproduce the experimental data of the neutron mul- 
tiplicity from the spherical region, that is, to reproduce the peak which is 
located at u n —8, where a n = A/12.0 is used. The neutron multiplicity in this 
case is shown by the dashed lines in Fig. 9. The black and gray lines denote 
the neutron multiplicity from the QF and FF processes, respectively. We can 
see that the value at v n = 6 is rather large, which is not consistent with the 
present experimental data. 

According to reference flfih . a n depends on the deformation. As we did in the 
previous section, we use a different value inside the fusion box and outside of 
it, a s v h = A/12.0 and a d ^ = A/10.0, respectively. We also plot this result in 
Fig. 9 as the solid lines. Due to the larger value of in the deformation space 
compared with a^ h , the emission rate of neutrons on the QF path decreases. 
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Fig. 10. The neutron binding energy B n of 266 Ds on the deformation space, (a) z — a 
(5 = 0) plane, (b) z — 6 (a = 0) plane. 



3.2-4 Neutron binding energy 5, 



Another parameter that strongly influences the decay width for neutron emis- 
sion is the neutron binding energy, B n , of the composite system. As shown 
in Eq. (7), B n influences the lifetime of neutron emission in the same manner 
as does a n . B n depends on the deformation of the nucleus. Figure 10(a) and 
(b) shows B n of 266 Ds on the z — a (6 = 0) plane and z — 5 (a = 0) plane, 
respectively. Here, we use the two-center shell model code f)l7l :ll8h. We can see 
that in the whole deformation space, the neutron binding energy fluctuates by 
about ±0.3 MeV comparing with that for spherical compound nucleus. 

We investigate the B n dependence of the neutron multiplicity. In the same 
way as in the case of a n , we use a different value of B n inside the fusion box 
{B^ jh ) and outside of it {B^). B^ h means the neutron binding energy of the 
compound nucleus. Figure 11 shows the neutron multiplicity for three different 
sets of B*>f t that is to say, = B s ? h } = B%> h - 0.3 and B^ = B^ h + 
0.3, as denoted by the solid, dashed and dashed-dotted lines, respectively. The 
black and gray lines denote the QF and FF processes, respectively. The larger 
value of B n causes the emission rate of neutrons to decrease. As a result, the 
distribution of the neutron multiplicity shifts to smaller multiplicity (to the 
left in Fig. 11). The large value of B n contributes the neutron multiplicity in 
the same manner as the large value of a n . 
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Fig. 11. The neutron multiplicity calculated using B*>f = B%> h ,B%> h - 0.3, and 
B^ ph + 0.3, as denoted by the solid, dashed and dashed-dotted lines, respectively. 
B^> h means the neutron binding energy of the compound nucleus. The black and 
gray lines denote the QF and FF processes, respectively. 

4 Summary 



We introduced the effect of pre-scission neutron emission into the three-dimensional 
Langevin calculation, that is to say, we combined the Langevin code with the 
statistical code. We applied our model to the investigation of the whole fusion- 
fission dynamical path in the reactions 58 Ni+ 208 Pb at £'*=185.9 MeV (jl). 

Our model is successfully applied to distinguish between the FF process and 
the QF process in terms of the pre-scission neutron multiplicity. We can see 
clearly the two components of the neutron multiplicity, which come from the 
QF process and the FF process. It is found that the peak for the large neutron 
multiplicity originates from the FF path, and on the other hand the peak for 
small neutron multiplicity comes from the QF path. The pre-scission neutron 
multiplicity has a strong correlation with dynamical path and is a powerful 
tool for investigating fusion-fission mechanism. 

We also discussed the sensitivity of the parameters of our model and showed 
clearly the role of each parameter. As the essential factors that influence the 
neutron multiplicity, we investigated the friction force 7, dissipation energy 
Ediss, level density parameter a n and neutron binding energy B n . 

This study will be the cornerstone for further study which develops the meth- 
ods of the classification between the FF process and the DQF process. The 
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DeMoN group already has measured the pre-scission neutron multiplicity in 



reactions, 48 Ca+ 208 Pb, 48 Ca+ 244 Pu and 58 Fe+ 248 Cm at E* ~ 40 MeV fl20j;|2l|). 
In the next study, we would like to analyze these reactions. 

The authors are grateful to Professor Yu. Ts. Oganessian, Professor M.G. Itkis, 
Professor V.I. Zagrebaev for their helpful suggestions and valuable discussion 
throughout the present work. The authors thank Dr. S. Yamaji and his col- 
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two-center parameterization. This work has been in part supported by INTAS 
projects 03-01-6417. 
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